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ABSTRACT 

We present an analysis of the 101 ks, 2007 Suzaku spectrum of the LINER 
galaxy NGC 1052. The 0.5 — lOkeV continuum is well-modeled by a power-law 
modified by Galactic and intrinsic absorption, and it exhibits a soft, thermal 
emission component below 1 keV. Both a narrow core and a broader component 
of Fe Ka emission centered at 6.4 keV are robustly detected. While the narrow 
line is consistent with an origin in material distant from the black hole, the broad 
line is best fit empirically by a model that describes fluorescent emission from the 
inner accretion disk around a rapidly rotating black hole. We find no evidence 
in this observation for Comptonized reflection of the hard X-ray source by the 
disk above lOkeV, however, which casts doubt on the hypothesis that the broad 
iron line originates in the inner regions of a standard accretion disk. We explore 
other possible scenarios for producing this spectral feature and conclude that the 
high equivalent width (EW ~ 185 keV) and full-width-half-maximum velocity 
of the broad iron line (v > 0.37c) necessitate an origin within d ~ 8r g of the 
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hard X-ray source. Based on the confirmed presence of a strong radio jet in this 
galaxy nucleus, the broad iron line may be produced in dense plasma near the 
base of the jet, implying that emission mechanisms in the centralmost portions 
of active galactic nuclei are more complex than previously thought. 

Subject headings: galaxies:individual - accretion, accretion disks - black holes - 
galaxies:nuclei - X-rays:spectra 



Introduction 



Since their first identification as an interesting class at optical wavelengths (IHeckman 



19801 ). the nature of Low Ionization Nuclear Emission Regions (LINERs) has remained elu- 
sive. Many such regions are associated with active galactic nuclei (AGN), though their X-ray 
luminosities are highly sub-Eddington and significantly below those of their Seyfert galaxy 
counterparts. The reason for this modest nuclear energy output remains unknown: possi- 
bilities include heavy absorption of the nuclear emission, comparably inefficient radiation of 
the accretion flow, or some combination of both. 

The notion of an advection-dominated accretion flow (ADAF) may apply in LINERs, 
wherein the ions lose ther mal contact with the electrons and onl y a small fraction of the dis- 
sipated energy is radiated (IRees et al.lll982l ; lNarayan &: Yilll994j ). Much controversy remains 
as to whether this picture is c orrect, however, and the debate has become more pointed since 
the report by lHo et al.l (119971 ) that some 20 — 30% of all galaxies are members of the LINER 
class. 

X-ray spectroscopy can help in determining whether the accretion disk is in an ADAF 
state or not. A LINER galaxy harboring an AGN produces a characteristic hard X-ray 
continuum. If the accretion flow close to the black hole conforms to the conventional thin 



disk archetype (jShakura fe Syunyaevl Il973l ). the cool gas orbiting in the inner disk should 
intercept and reprocess a significant fraction of the continuum, also producing the Fe Ka line 
(among many fluorescent emission lines from other species as well). This spectral feature 
should have a broadened, skewed profile created by the relativistic properties of the spacetime 
close to the black hole, and a significant Comptonized reflection component should be seen 
in the continuum above 10 keV as well. If, however, the accretion flow in this region takes the 
form of an ADAF, the ion temperature is so high that the iron ions should be totally stripped 
of electrons. Furthermore, the trapped energy will puff up the disk, reducing its density and 
optical depth. Both of these conditions render the reflected emission, including the emission 
of Fe Ka, virtually non-existent in this region. In this case, only contributions to the Fe Ka 
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line from other locations (e.g., from a disk wind, the outer disk, broad line region or putative 
molecular torus) will manifest, and the line may likely have a much narrower profile as a 
result. 

The LINER galaxy NGC 1052 is an interesting case in which to study the dynamics of 
the accretio n flow. The AGN is housed by a nearby elliptical host galaxy with a redshift of 
z = 0.0049 (IKnapp et al.lll978l ). implying a distance of 20.7 M pc using WMAP co smology. 
The source has a well-studied LINER optical spectrum (e.g., iGabel et al.l (120001 )). and is 
also reasonably bright in the X-ray band for a LINER, with an a yerage 2 — lOkeV flux 
of F2-10 ~ 5 x 10~ 12 erg cm -2 s _1 . Using the M — a relat io n of iTremaine et al.l (120021 ) 



and the stellar velocity dispersion of iNelson fe Whittle! (119951 ). IWoo fe Urryl (120021 ) derived 
a black hole mass of Mrh = 1-54 x 1O 8 M for NGC 1052, and a bolometric luminosity 
of Lboi = 6.92 x 10 43 ergs _1 using flux integration over its spectral energy distribution. 
Combining these results with the derived Eddington luminosity yields an estimated accretion 
efficiency for the black hole in NGC 1052 of 77 < 0.004, in Eddington units. This is consistent 
with general results from LINERs. The so urce is known to be heavily absorbed based on the 
presence of H2O megamasers in its core ( 
of a radio jet also present in the galaxy ( 



Claussen et al 



Bra atz et al.l 119941) . which lie along the direction 



19981 ). Recent VLBI studies have 



constr ained the jet inclination angle to lie between 57 — 72° to our line of sight (IKadler et al. 



2004bl ). implying that the jet emission is de-beamed, and suggesting that the bulk of the 
X-ray emission is not produced by the jet. 

The observed X-ray spectrum of NGC 1052 is quite flat (for the power-law component 
of the emission T ~ 0.2 — 1.7, depending on the data and models used ) and appears to be 



dominated at lower energies by thermal an d /or photoionized emission (jWeaver et al.l 11999 



Guainazzi et al.l |2000| ; IKadler et al.l l2004al ) . The source of this soft emission is likely the 
interaction of the jet with the surrounding ISM (IKadler et al.l l2004al ) . while the power-law 
component at harder energies is thought to originate from nuclear emission. Radio continuum 
and RXTE monitoring both show that the source varies strongly on time scales of weeks to 
months. The 2 — lOkeV flux, in particular, varies between ~ 4 - 9 x 10~ 12 erg cm -2 s^ 1 , 
and though there are large uncertainties in the RXTE data, the spectrum appears to switch 
from softer to harder states, corresponding to changes in overall source flux (Kadler et al. , 
in prep.). 



A prominent Fe Kg line has been pre viously observed in the source by ASCA (jWeaver et al 



1999), BeppoSAX (iGuainazzi et al 



dra (IKadler et al.ll2004bl ; iRos et al. 



2000), XMM-Newton (Kadler et al. , in prep.), and Chan- 



20071 ). Of these observations, a broadened component to 



the Fe Ka line was seen with BeppoSAX and XMM-Newton. Given the low a ccretion rate of 
NGC 1052 — a scenario which has long been thought to produce ADAFs, e.g.. lNarayan fe Yi 
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( 119941 ) — it is natural to speculate on the origin of the broad iron line. The spectral state 
of NGC 1052 is in many ways analogous t o the low/hard sta te observed in X-ray binaries 
(XRBs). Contrary to the ADAF paradigm. Miller et al.l (120061 ) reported that XMM-Newton 
observations of the relativistically broad Fe Ka line in the low/hard state of XRB GX 339- 
4 showed a standard thin disk remaining at or near to the innermost stable circular orbit 
(ISCO), at least in bright phases of the low/hard state. The authors therefore argued that 
potential links between the inner disk radius and the onset of a steady compact jet, and the 
paradigm of a radially recesse d disk in the low/hard state, do not hold universally. Note, 
however, that iGierlihski et al.l (120081 ) cautioned against using only the shape of the Fe Ka 
line to determine the radius of the ISCO, suggesting that these results can be significantly 
underestimated unless irradiation from the energetically dominant hot plasma (along with 
other spectral effects) is taken into account. If the broad iron line in NGC 1052 does origi- 
nate from the accretion disk, however, then this s pectral feature, coupled with detailed VLBI 
observations of the inner jet (IKadler et al.ll2004al Jbl). may offer a unique opportunity to study 
the connection between the accretion flow and jet production in AGN. 

A li nk between accretion and jet activity in AGN has been no ted in the radi o galaxies 
3C120 fcataoka et al.l 120071 : iMarscherl lioOft iMarscher et al.l l2002h and 3C 111 jMarscher 
20061 ). While Marscher et al. have noted dips in the X-ray light curves of these sources 
correlated with ejections of superluminal knots in the radio jet, suggesting that a portion 
of the inner disk is disturbed during an ejection event, these results fall short of mapping 
out the physical structure of the disk at the time. However, Kataoka et al. observed a broad 
Fe Ka line in 3C 120 with Suzaku in both the high and low flux states of the source. This 
discovery demonstrates the viability of studying the accretion flow and the jet activity of 
radio-loud AGN simultaneously, and will hopefully pave the way for new insights into the 
fundamental link between these two processes. 

Here we discuss the analysis of our 2007 Suzaku observation of NGC 1052, which offers 
the highest quality X-ray spectrum of this galaxy to date and gives us an unprecedented, 
simultaneous view of the 0.5 — 35keV energy range. NGC 1052 is one of a small handful of 
AGN sufficiently bright in both the X-ray and radio bands for tests of a jet/disk connection. 
Additionally, it is the only source for which we will have a full package of multi-epoch X- 
ray spectral data, X-ray and radio light curves, and VLBA imaging within a common time 
window and with an analysis by a single collaborative group. A forthcoming paper will 
provide full details of this ongoing campaign, placing our Suzaku results in the context of 
prior X-ray and ongoing radio analysis of this source (Kadler et al. , in prep.). 

Our Suzaku observation of NGC 1052 and the data reduction will be detailed in Section 
2. We discuss variability over the course of the observation in §3. Our approach to the 
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spectral fitting of the time-averaged data and the results of this analysis for both the XIS 
and the HXD/PIN instruments are presented in §4. We examine the Fe Ka region at length 
in §4 as well. A comparison with results from previous epochs appear in §5, and we discuss 
our results in context in §6. Our conclusions comprise §7. 

2. Observations and Data Reduction 

We observed NGC 1052 for a total of ~ 101 ks from 16-18 July 2007 with Suzaku in 
the HXD nominal pointing position. Reduction of the Suzaku data followed the procedures 
outlined in §4.7 and §4.8 of the Suzaku ABC Guide, available online^. Calibration files used 
were the latest version as of the time of this writing (January 2009). This calibration was 
incorporated into the X-ray Imaging Spectrometer (XIS) data reduction using the xispi 
task on the unfiltered data. Once cleaned, the event files for the XIS 0, XIS 1 (XIS 2 became 
defunct as of Nov. 2006) and XIS 3 detectors were loaded into xselect for reprocessing. 
Images were extracted that enabled identification of source and background regions. As per 
the recommendations in the ABC Guide, we employed a circular source region ~ 250" in 
radius, and maximized our background region area while avoiding the calibration sources 
at the corners of each detector chip. After spatial filtering was applied, spectra and light 
curves were extracted for each of the operational three detectors. Response matrices and 
ancillary response files were generated using the xisresp task, and the spectra, backgrounds 
and responses for the front-illuminated chips (XIS and XIS 3) were co-added to increase 
the signal-to-noise using addascaspec. The spectra, backgrounds and responses were then 
rebinned by a factor of eight (to 512 channels) using the rbnpha and rbnrmf tasks to speed 
spectral fitting. Finally, the spectra, backgrounds and responses were grouped together using 
the grppha task. Only spectral channels with with a minimum of 20 cts/bin were used for 
fitting in order to ensure the validity of the chi-squared statistic. This corresponded to an 
energy range of 0.7 — 9keV for the XIS 0+3 data (32113 total counts) and 0.5 — 7keV for 
the XIS 1 data (16146 total counts). 

Suzaku did observe NGC 1052 with the silicon diode PIN instrument of the Hard X-ray 
Detector (HXD), but the GSO crystal scintillator instrument did not have enough counts 
to provide useful data. The PIN data were reduced using the "tuned" calibration files for 
the non-X-ray background (NXB) and the response files that had been generated for the 
period surrounding the observation. These files were downloaded from the Suzaku Guest 
Observer Facility (GOF) for use. We first needed to ensure that only the portion of the 



http: / /heasarc / docs / suzaku / analysis /abc/ nodel .html 
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NXB coincident with our observation was used in processing, so we merged the good-time 
interval (GTI) of the NXB with that of our screened event file to yield a common GTI using 
mgtime. With xselect, we then filtered the data and background using this common GTI 
and extracted our spectra and light curves for the source and NXB. We corrected for dead 
time in the observation using the hxddtcor task, and increased the exposure time of the 
NXB data by a factor of ten to compensate for the inflated event rate, which had been 
similarly multiplied in an effort to suppress Poisson noise. As for the XIS data, only spectral 
channels with with a minimum of 20 cts/bin were used for fitting, which limited our energy 
range to 12 — 35 keV. After reduction and filtering, the PIN spectrum had 35705 total counts 
over this energy range. 

The cosmic X-ray background (CXB) has not been taken into account in PI N obser 



vations, so to model it appropriately, we simulated the CXB spectrum in xspec f lArnaud 



19961 ) using a power-law model of T = 1.29 with a high energy cutoff and foldE= 40keV 
(initial parameters are given in the ABC Guide). We then folded this model through the 
response file for the PIN flat field for an exposure time of 10 6 s to generate the simulated 
CXB data. This spectrum was then read back into xspec with the PIN response for our 
observation and the original model was refitted over 12 — 35 keV to yield the actual param- 
eters of the CXB during this time. In our case, the simulated CXB spectrum contributed 
a count rate of ~ 2.28 ± 0.02 x 10 _2 ctss _1 (~ 5%) to the total X-ray background from 
12 — 35keV. In modeling the PIN data with xspec, we applied the best-fit CXB parameters 
as a constant with each fit. The CXB followed a power-law form with a high-energy cutoff: 
T = 1.30 ±0.19, K po = 8.13 ±3.07 x 10~ 4 phcm- 2 s" 1 , foldE = 40.53 ± 15.93 keV. The PIN 
instrument detected the source above the (NXB±CXB) background at a > 10cr level in our 
observation. 



3. Timing Analysis and Variability 

No significant changes were seen in the count rates of any of the XIS instruments or 
the PIN instrument over the course of the observation. Fig. [1] depicts the light curves from 
XIS 0, XIS 1 and XIS 3, as well as the combined XIS light curve vs. that of the PIN. While 
small variations in count rate (over ~ 10 4 s timescales) were seen in both data sets, neither 
varied more than a factor of two from its baseline flux throughout the observation. This 
flux was roughly 0.07 ctss -1 for the PIN and 0.38 ctss -1 for the combined XIS data (0.16 for 
XIS 0, 0.14 for XIS 3 and 0.21 for XIS 1). Fitting the combined XIS light curve to a constant 
model resulted in x 2 l v — 227 7236 (0.96) (where v denotes the number of degrees of freedom 
of the model: v = # detector channels - # free model parameters), and an uncertainty of 
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< 0.01 ctss" 1 on the 0.38 ctss" 1 count rate. The PIN light curve fit to a constant was not 
as good a fit at x 2 / u = 351/245 (1.43), with the uncertainty in its 0.07 ctss -1 count rate at 

< 0.005 ctss -1 . We also computed the hardness ratio vs. time for the combined XIS data, 
which is consistent with unity as one might expect for a source without significant variability 
over the course of the observation: HR = 1.02±0.02, with x 2 / u = 101/71 (1.42). Examining 
the hardness ratio vs. soft band count rate, we note a significant anti-correlation between 
the two: NGC 1052 becomes harder as the source flux decreases. This variation is best fit 
with a linear model: HR = m(flux) + b, where HR denotes the hardness ratio, "flux" is the 
soft band count rate, m is a multiplicative term and b is a constant term. For the combined 
XIS data, m = -8.35 ±2.31 and b = 1.92 ±0.25, with X 2 /v = 66/71 (0.93). We have plotted 
the hardness ratio vs. time and vs. soft count rate in Fig. [2j 

This hardening of the source with decreased count rate demonstrates in a model- 
independent way that the spectral variability of of NGC 1052 in X-rays is consistent wit h 



the behavior observed in ot her Seyfert or Seyfert-like AGN (e.g., iPapadakis et al.l (120021 ): 



Markowitz fe Edelsonl (120041)'). and contrasts with the behavior observed in jet-dominated 



AGN (e.g., IZhang et al.l ((2006)). In other words, the anti-correlation of the hardness ratio 
with the source flux confirms that the bulk of the X-ray emission in NGC 1052 is not pro- 
duced by the jet. We have plotted the hardness ratio vs. time and vs. soft count rate in 
Fig.H 

The overall 2 — lOkeV source flux did not vary significantly, averaging to i^-io = 
5.37 x 10 -12 erg cm -2 s -1 (absorbed). This was split between the continuum and the Fe 
K lines, where the continuum flux averaged F cont = 4.86 x 10 -12 ergcm -2 s -1 and that of 
the combined broad and narrow Fe Ka lines averaged Fp e = 5.10 x 10~ 13 erg cm -2 s -1 . 
The average unabsorbed 2 — lOkeV flux at the distance of NGC 1052, F 2 -io = 8.96 x 
10 -12 erg cm -2 s -1 , works out to an intrinsic luminosity of L2-10 — 4.60 x 10 41 ergs _1 using 
cosmological parameters consistent with the WMAP5 ACDM-model results, within errors: 
H = 70kms -1 Mpc -1 , q = 0.00 and A = 0.73. These results, and the models used to 
derive them, are presented in more detail in §4. 



4. Spectral Analysis 

We began our analysis of the 2007 Suzaku spectrum of NGC 1052 by examining the data 
from the operational XIS detectors. We use the good signal-to-noise 0.5 — 9.0 keV data to 
extract information about the underlying continuum, any soft excess emission, the Fe K line 
complex and any intrinsic absorption in the system. NGC 1052 is thought to exhibit evidence 
for an AGN disk (the broad Fe Ka line) and jets (radio observations), as well as substantial, 
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NGC 1052 Bin time: 5000. 

Bin lime: 5000. s 




Fig. 1. — Left: Light curves from the three operational XIS instruments for the 2007 Suzaku 
observation of NGC 1052. Right: Light curve of the combined XIS data (black) as well as that of 
the PIN (red) over the course of the observation. All light curves shown are background-subtracted, 
including CXB subtraction for the PIN data. Note the overall lack of significant variability of the 
source. 




Fig. 2. — Left: Hardness ratio plotted against time; Right: hardness ratio plotted against the soft 
flux over the course of the 2007 Suzaku observation of NGC 1052. Though there is no significant 
change in hardness ratio vs. time during the observation, note that as the soft flux drops slightly, 
the source becomes harder. 



- 9- 



extended soft emission. As such, we expected the spectrum to be a conglomerate of all of 
these features since Suzaku lacks the spatial resolution necessary to physically separate them. 

We also analyzed the HXD/PIN spectrum from 12 — 35keV in order to learn more 
about the nature of the continuum and to constrain the amount of reflection seen in the 
source. Because a bro ad iron line component has been observed in NGC 1052 in the past 
(iGuainazzi et al. I hoool ). we expected that residual emission above a power-law should remain 
at high energies due to the presence of the so-called "Compt on hump." This fea t ure h as 
been noted in many AGN with broad iron lines (for a review, see lReynolds fc Nowakl (120031 )). 
representing Compton down-scattering of the hard X-ray photons by the material in the disk. 
This Compton hump typically peaks at ~ 20 — 30 keV, and is thought to go hand-in-hand 
with the presence of a broad iron line, as both spectral features arise from the same physical 
process of reflection. 



4.1. XIS Continuum 



Before the iron line region could be considered, it was necessary to first accurately model 
the underlying continuum of NGC 1052. We began by ignoring the energy range where the 
iron line is thought to be important (3 — 7keV) and modeled the 0.5 — 3 and 7 — 9keV 
energy bands wit h a simple power-law modified by Galactic photoabsorption (Ah = 2.83 x 



10 20 cm- 2 



as per 



Kalberla et all ((2005J)). The result was a poor fit: ^jv = 840/361 (2.33). 



Large residuals remained at both high and low energies (see Fig. [3]), and a very hard photon 
index for the power-law component was seen: T = 0.34 ± 0.12. To account for differences 
in calibration between the XIS 0+3 and the XIS 1 detectors, we multiplied each data set 
by a constant cross-normalization factor. This constant was held at 1.00 for XIS 0+3, and 
achieved a best-fit value of 0.94 ± 0.01 for the XIS 1 data. 

An apparent soft excess was clearly visible below ~ 1.5 keV ( see Fig. [3]). To m itigate 
this residual feature we incorporated a thermal mekal component (IMewe et al.lll985l ). The 
fit improved dramatically to ^jv = 381/358 (1.06) (Ax 2 /Av = -459/ - 3, or an improve- 
ment in the global goodness-of-fit by 459 chi-squared points for an additional 3 degrees of 
freedom), with kT = 0.6 ± 0.05 keV and a solar abundance for the me kal component. This 
component has bee n noted in preyious X -ray analyses of NGC 1052 by lKadler et al.l (l2004bl ) 
with Chandra and I Weaver et al.l (119991 ) with ASCA (hereafter W99), and the temperature 
and flux obtained from our Suzaku fit are comparable to both values in both previous works. 
This thermal emission most likely origi nates in extended em ission around the nucleus from 
the interaction of the jets with the ISM (IKadler et al.ll2004bl ). Residuals still remained at low 
energies, however, so we added in a second mekal component (with an abundance equal to 
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Energy (keV) 

Fig. 3. — The XIS spectrum of NGC 1052 fit with a power-law modified by Galactic photoabsorp- 
tion over the 0.5 — 3 and 7 — 9keV bands uninfluenced by the fluorescent Fe Ko line. The fit is 
poor: x 2 l v = 840/361 (2.33). Note especially the large residuals at both the hard and soft ends. 
Black crosses represent the combined XIS 0+3 spectrum while the black line represents the model. 
Red crosses and line represent the XIS 1 data and model. The lower panel of the figure shows the 
data-to-model residuals, with the green line representing a perfect fit with a data-to-model ratio of 
one. 
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that of the first) to try to account for these features. Although marginal improvement in the 
overall goodness-of-fit was seen (Ax 2 /Ai/ = —12/ — 2), the addition of this second thermal 
component rendered the parameters of both mekal components unconstrained in an error 
analysis. We therefore removed the second thermal component. Though modest residuals 
remained at select energies, no statistically significant improvement in fit was achieved with 
the addition of discrete Gaussian components, as one might expect in a photoionized plasma, 
or with allowing the abundance of the mekal component (s) to vary. Alternatively, we also 
tried to account for this excess soft emission with a second power-law, with and without 
further absorption. Both models exhibited significantly poorer fits than that of the single 
mekal component. 

Because NGC 1052 is also classified as a Seyfert 2 galaxy, it is reasonable to assume 
that some intrinsic absorption exists in or around the nuclear region. Such absorption could 
play a significant role in the spectral curvature seen, so it must be properly accounted for in 
the continuum model. We employed a partial-covering absorber (zpcf abs) intrinsic to the 
source to model this component, applying it only to the power- law and other components of 
nuclear AGN emission and not to the extended mekal components. A statistically significant 
improvement in fit was seen: x 2 jv = 335/356 (0.94), with Nn = 1.08 ±0.15 x 10 23 cm -2 and 
a covering fraction of 84 ± 2%. This marked an improvement in fit of A% 2 /Az/ = —46/ — 2. 
With the inclusion of this intrinsic absorption, the power-law photon index increased to 
T = 0.80 ± 0.09. This was a significant increase in photon index from the fit without 
intrinsic absorption. 

Physically, this continuum model and its parameters were roughly comparable to those 
which were fit f or the 1996 ASCA obse rvation of the source (W99) and the 2000 BeppoSAX 



observation by iGuainazzi et al.l (120001 ) — hereafter GOO — although the photon index of 
the power-law is considerably harder and the intrinsic absorption less dense in this 2007 
observation. A model with a dual neutral absorber modifying intrinsic plus scattered power- 
law components as in W99 and GOO did not improve the fit, statistically; moreover, the 
parameter values of the second power-law and absorber could not be constrained. 

Even after modeling the neutral absorption in the system, some residual spectral cur- 
vature still remained around 2 — 3 keV, indicating that some further component of absorp- 
tion likely remained unmodeled. We applied the simple ionized absorber model absori 



(IZdziarski et al.lll995l ) to the nuclear emission and noted a marked visual improvement to 
the fit from 2 — 3keV, as well as a small improvement in the global goodness-of-fit to 
X 2 fv = 325/354 (0.92), an improvement of A% 2 /Az/ = —10/ — 2. This absorber was reason- 
ably thick and moderately ionized, with a column density of Nr = 1.37 ± 0.02 x 10 22 cm -2 
and an ionization parameter of £ = 68±7ergcm _1 s" 1 . With the addition of this component, 
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the power-law photon index increased to T — 1.50 ± 0.02. This value lies within reasonable, 
expected physical limits for the X-ray continuum of an AGN. Our base continuum model 
therefore included a power-law modified by intrinsic absorption from both a cold, patchy 
absorber and a moderately ionized absorber of smaller column density, along with extended 
thermal emission. All model components were affected by Galactic photoabsorption. 



4.2. The Fe Ka Line 

Having fit the continuum, we then took the energies from 3 — 7keV back into consid- 
eration and attempted to model the prominent residuals remaining, highlighted by a strong 
emission feature centered at 6.4 keV in the rest frame that was assumed to be neutral Fe 
Ka (see Fig. HJ top). We began by holding the continuum parameters constant, except for 
normalizations and the power-law photon index, and refitting. The strong residual at 6.4 keV 
remained, appearing to have a broadened red wing associated with it that extended down 
to ~ 4keV. For this fit, x 2 /v = 711/631 (1.13). 

To account for the strong emission feature of 6.4 keV, we initially assumed some core 
contribution from neutral Fe Ka outside of the inner accretion disk in distant material. This 
would produce a relatively narrow emission line, being outside the region where relativistic 
broadening becomes important. Inserting such a line (a = 5eV, intentionally less than the 
resolution of the XIS in order to make the line truly narrow) into the model yielded a vast 
improvement in fit to x 2 /" = 584/629 (0.93) for Ax 2 /Av = -127/ - 3 over the case with 
no line present. The energy of the line was constrained in the fit to E — 6.39 — 6.41 keV, 
indicative of neutral iron, and the feature had an equivalent width EW = 88 — 134 eV. These 
parameter ranges are quoted at 90% confidence for one interesting parameter; this confidence 
interval applies to all parameter ranges quoted in this Section. The broad residuals centered 
on 6.4 keV still remained, however (see Fig. HI middle). 

We attempted to model these broad residuals three separate ways: with a broad Gaus- 
sian, with a diskline compone nt representing emission from the inner disk around a Schwarzschild 



black hole (IFabian et al.lll989l ). and with a laor compo nent repres enting emission from the 



inner disk around a maximally-spinning Kerr black hole (lLaorlll99ll ). Each line was centered 
at 6.4 keV in the rest frame, corresponding to neutral Fe Ka. The Gaussian component 
improved the fit by A% 2 /Az/ = —13/ — 2 to x 2 l v = 571/627(0.91), and was quite broad: 
a = 642 — 852 eV. Its equivalent width of -EWbroad = 166 — 236 eV was slightly smaller than 
that of previous observations, e.g., W99 and GOO, in which EWi, Toa d ~ 300 eV. Relaxing the 
energy constraint of the line resulted in E = 5.41 — 5.88 keV, with a = 780 — 1300 eV and 
EWfroad — 359 — 390 eV, though significant residuals still remained on the red wing of the 
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Fig. 4. — Top: The residual emission feature remaining after fitting the continuum of the NGC 1052 
Suzaku/XIS data, as per §4.11 Note the prominent peak centered at 6.4 keV in the rest frame 
and the excess emission extending down to ~ 4keV. x 2 l v = 711/631 (1.13). Middle: Residuals 
remaining after fitting only the narrow Fe Ka core. Evidence fo r a broader feature is clear, x 2 l v = 
584/629(0.93). Bottom: Residuals remaining after a laor line (|Laorlll99ll ) was used to model the 
broad Fe Ka emission line as well, x 2 l v = 561/625 (0.90). 
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line. We then turned to diskline and laor in the hope of modeling the underlying broad 
feature, recognizing that these models assume that the broad emission originates in the inner 
accretion disk, which may not be the case if NGC 1052 harbors an ADAF. 

The diskline model yielded x 2 / u = 574/625 (0.92), which is not significantly different 
from the Gaussian fit (A% 2 /Az/ = +3/ — 2), statistically, though it did account for more 
of the residuals below 6 keV. VLBI observations of the radio jet on sub-parsec scales have 
constrained t he inclination angle of the accretion disk to be between 57 — 72° relative to our 
line of sight (IKadler et al.ll2004bl ). so we applied these constraints to the model. Our best 
fit yielded an accretion disk emissivity index of a = 1.69 — 3.69 (where the disk emissivity 
e oc r~ a ) and an inner disk radius of r- m < 26 r g (where r g = GM/c 2 ), but the inclination 
angle could not be further constrained. We held the outer radius constant at r out = 400 r g , 
since it was unable to be constrained and our emissivity index dictated that this radius safely 
encompasses all of the relevant emission from the disk. Fixing the inclination angle to an 
average value of 65° and the emissivity index of the disk to a "typical" Seyfert AGN value 
of a = 3 yielded no significant improvement in fit or change in the other model parameter 
values, and only resulted in a constraint on the inner radius of disk emission of rj n = 1 1 —52 r g . 

The laor model provided the best statistical fit at \ 2 l v — 561/625 (0.90), an improve- 
ment of Ax 2 / Azy = -13/0 over the diskline model and A% 2 /Az/ = —10/ — 2 over the 
broad Gaussian line model. The laor model also appeared to be most effective at reducing 
the residuals between 4 — 7keV (see Fig. HI bottom). Keeping the same disk inclination 
constraints mentioned above for the diskline model, we again could not further constrain 
this parameter, but we noted that the emissivity index of the disk was fairly centrally con- 
centrated at a = 1.53 — 5.10, while the inner radius of emission is only mildly constrained 
within r in = 8.93 — 45.10 r g . The equivalent width of the broad iron line here was found 
to lie within .EWbroad = HO — 297 eV. Interestingly, when we relaxed the constraint on the 
inclination angle of the disk we obtained a best-fit value of i = 45 ± 5°; this is roughly a 
~ 2cr deviation from the radio constraints. Fixing the inclination angle and emissivity index 
as above for the diskline case worsened the global goodness-of-fit by 22 chi-squared points 
for 2 fewer degrees of freedom (A^ 2 /Az/ = +22/ + 2) and no significant improvement on the 
constraint for the inner radius of disk emission: r in = 13 — 44 r g . 

The laor fit was a slight improvement over both the Gaussian and diskline models, 
though it should be noted that the uncertainty of the r- m parameter effectively renders it 
impossible to meaningfully distinguish between the diskline and laor models based on 
the Fe Ka profile alone. If the laor model fit could better constrain the inner disk radius 
to r in < 6r g (indicating that an optically-thick disk extends down to small radii), it would 
suggest that the black hole at the center of NGC 1052 is rapidly rotating, as many theories 
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expec t for a source that p o wers r elativistic jets (e.g.. lBlandford fc Znajekl (119771 ) ; lMeier et al. 
( 120011 ); lHawley fc Krolikl (120061 )). We cannot make this statement unequivocally based on 
our data, however. The data fit with this laor model are shown in Fig. The final best-fit 
model parameters and their errors, incorporating the HXD /PIN data as well, are presented 
in the following Section within Tabled] as Model 1. 

Absorption has often been suggested as an alternative to broad iron e mission when 
attem pting to model enhanced spectral curvature in the 4 — 7keV range (e.g.. iKinkhabwala 
( 120031 )). To address this possibility, we compared our best laor fit above to a model in which 
we removed the broad iron line and replaced it with a second absorbed power-law component. 
The resulting fit yielded \ 2 > = 582/627(0.93), as opposed to the ^jv = 561/625(0.90) 
achieved with the final laor model (A% 2 /Az/ = +22/ + 2), and left significant residuals 
around 6.4 keV and below. Moreover, the parameter values for both the partial-covering 
absorber and the power-law components are not constrained in the fit. The models employing 
broad emission lines unquestionably provide better statistical and visual fits to the data. 



4.3. Addition of the HXD/PIN Data 

Adding the HXD/PIN data on to the combined XIS 0+3 and XIS 1 spectra, we noticed 
that the signal-to-noise ratio deteriorated significantly above ~ 35keV and below ~ 12keV. 
As such, we restricted our energy range for this detector to 12 — 35keV, applying a constant 
cross-normalization factor as we did with the XIS 1 data set. In this case, we used a 
parameter value of 1.18 in keeping with the value expected for the HXD nominal pointing 
position. 

We considered the combined data refitted with Model 1 (photoabsorbed power-law, 
a partly-ionized absorber and a partial-covering neutral absorber, soft thermal emission, 
narrow and broad Fe Ka components), including no reflection component above 9 keV. Fig. [6] 
shows the data, model and residuals from this fit, which yields y 2 jv = 660/685 (0.96). The 
inclusion of the high-energy data left the 2 — 9 keV parameters effectively unchanged from 
their values in the fit without the PIN data. If a significant flux contribution from reflection 
is present, as one might assume if the broad iron line has an origin in a relativistic accretion 
disk, we would expect to see some excess in the emission over 9keV relative to Model 1. No 
such feature was detected, however. 

Table [1] shows the full 0.5 — 35keV fit with 90% confidence errors. 
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Fig. 5. — The best-fitting model for the XIS data. This model includes a power-law continuum 
modified by intrinsic absorption from a partial-covering neutral material as well as a mildly ionized 
gas, extended thermal emission, and both broad and narrow components of neutral Fe Ka. The 
spectrum is also modified by Galactic photoabsorption. x 2 l v = 561/625 (0.90). The red, black and 
green colors/lines are the same as in Fig. 1. 
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Fig. 6. — HXD/PIN data included from 12 — 35keV and refitted with Model 1, which does not 
include any type of continuum disk reflection. Inclusion of this data does not significantly alter the 
model parameter values from the 0.5 — 9keV fit. Though a broad line is still required at 6.4 keV, 
no accompanying evidence of reflection from the disk above 9 keV is seen. Red and black colors 
follow the designations of Fig. 1. In the top panel, the green crosses and line represent the PIN 
data and model, respectively. In the bottom panel the horizontal green line still represents a data- 
to-model ratio of one, while the green crosses show the residuals between the PIN data and model. 
x 2 / v = 660/685 (0.96). 
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4.4. Alternative Continuum Models 

Though no evidence for Comptonized reflection is apparent to the eye above 9 keV, it is 
nonetheless important to try to place upper limits on its possible contribution to the overall 
spectrum. With this in mind, we utilized two popular, well-regarded models of reflection 
from an accretion disk around a black hole in an attempt to quantify the amount of reflection 
that might be present in the spectrum of NGC 1052. 



The pexrav (or pexriv) model of iMagdziarz fe Zdziarskil (119951 ) is one of the most 



widely used public models for AGN emission which incorporates reflected continuum from 
a neutral (or partly ionized) accretion disk. While it represents a more physical model than 
the simple power-law alone, it does not include line emission from the disk (e.g., Fe Ka) and 
it is difficult, in practice, to constrain the three most important parameters simultaneously: 
photon index, reflection fraction and disk inclination angle. 

Adding a pexrav component to our power-law continuum (and keeping the inclination 
angle of the disk constant at the laor value, the power-law cutoff at 100 keV, and the 
abundances fixed at their solar values), we found x 2 / v — 660/683 (0.97), only a very marginal 
change in fit from that of Model 1 over the 0.5 — 35keV energy range (A% 2 /Az/ = 0/ — 2). 
The photon index fit between T = 1.45 — 1.56 while the reflection fraction was constrained 
to be quite small: i? re fl < 0.01. As we expect, based on the residuals from fitting Model 1, 
reflection is not a significant contributor to the overall observed spectrum. Though a pexriv 
component was also fit, the ionization parameter was unable to be constrained, so we elected 
to assume a neutral disk. In this fit, the emissivity index of the laor line for Fe Ka had very 
similar constraints to those found in Model 1: a — 1.74 — 4.56. Also, to 90% confidence, the 
inner radius of emission in the disk was constrained to r in < 45 r g , as we found in Model 1. 
The equivalent width of the broad Fe Ka line in this model also remained consistent with 
that found in Model 1, at EW = 113 — 293 eV. The pexrav model will hereafter be known 
as Model 2. Its final parameter values and their errors are also presented in Table [H 



The ref lion model of iRoss fe Fabian! (120051 ) is another widely available code used to 



describe the reflected spectrum of an accretion disk around a black hole. The advantage 
of this model over the pexrav model is that it is self-consistent: fluorescent emission lines 
from many ionized species are included in addition to the continuum, most notably those 
of the Fe K complex. The ionization parameter, incident power-law photon index and iron 
abundance are specified by the user. The disadvantage is that photon indices below T = 1 
are outside the allowed parameter range of the model, and there is not a parameter for the 
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reflection fraction, as in pexrav. The user is therefore obliged to estimate it using, e.g., 



Here ifrefl an d K po denote the normalizations of the ref lion and power-law components 
from our best-fit spectral model, £ is the best-fit ionization parameter characterizing the 
disk reflection (where £ = 30 denotes a neutral disk in the model), and F refl and F po are the 
total fluxes contained in the ref lion and power-law components for unity normalization 
over the full wavelength range (0.001 — 1000 keV). 

We began by removing the laor line from the Model 1 fit and adding in a reflion 
component, with the iron abundance held at the solar value and the photon index linked 
to that of the power-law. As expected, this fit mimicked that of the continuum plus a 
narrow 6.4 keV Gaussian, leaving residuals on the low-energy end of the line. Statistically, 
X 2 / v = 679/687 (0.99). Allowing the iron abundance to fit freely yielded a small improvement 
in fit, with x 2 l v — 667/686 (0.97) for Fe/solar < 3.53, but doing so rendered the ionization 
parameter and normalization unconstrained. As such we elected to keep the iron abundance 
constant at the solar value. 

If NGC 1052 harbors a rapidly spinning black hole, as our best laor fit seems to indicate, 
we must take into account the relativistic effects on emission originating from the inner 
accretion disk. We incorporated this smearing of spectral features using the laor-based 
convolution model kdblur. Our best fit was x 2 l v = 660/684 (0.96), comparable to our best 
laor and pexrav fits (A% 2 /Az/ = 0/ — 1 and A% 2 /Az/ = 0/ + 1, respectively). However, we 
were unable to constrain the emissivity index of the accretion disk, though we did obtain 
constraints on the inner radius of emission that are consistent with those found in the laor 
model: r in < 45 r g . The estimated reflection fraction is quite low at i? re fl < 0.006, consistent 
with the results of the pexrav model fit and partly attributable to the moderately high 
ionization parameter (£ < 111) in the reflion model. The kdblur (ref lion) model will 
hereafter be known as Model 3. For comparison, the model components for the Models 1-3 
fits are plotted in Fig. and their best-fit parameters and errors are listed in Table [TJ 

Models 2-3 constrain the amount of reflection from the disk to be _R re fi < 0.01, consistent 
with the lack of a discernible Compton hump above 10 keV in the HXD/PIN data. Moreover, 
as we describe in detail later in SJ6], if the broad iron line component results from fluorescence 
in the disk caused by incident hard X-ray emission, then the equivalent width of the Fe Ka 
line should be directly related to the normalization of the reflection component in our model 
as compared to that of the power-law. Yet this is not the case in our data. Based on the low 
normalizations for reflection seen in both Model 2 and Model 3, we should expect a broad 
Fe Ka equivalent width of EW ~ 80 eV, much less than the observed EW ~ 185 eV in both 
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Models 1-2. In other words, our broad iron line does not appear to have a correspondingly 
strong amount of continuum reflection associated with it. It is possible that this "missing" 
reflection is actually present to some greater degree in the data, but is perhaps drowned out 
by the power-law continuum. It is also possible that reflection is simply absent in this source. 
This scenario presents several intriguing possibilities for physical interpretations, which we 
discuss in $61 



5. Comparison with Previous Results 

Here we consider our best 2007 Suzaku results (Model 1) against those obtained in 
the 1996 ASCA observation by W99 and the 2000 BeppoSAX obser vation by G00. We do 



not el ect to discuss individually either th e 1999 ROSAT ^result s from iGuainazzi fc Antonelli 



( 119991 ) nor the 2000 Chandra results from lKadler et al.l (j2004al ) due to the inferior signal-to 



noise of these observations. 

Neither W99 nor GOO presented a light curve for their data, but both authors failed 
to detect any significant variability over the course of their observations. W99 fit a time- 
averaged 2 — 10 keV flux of F2-10 ~ 8 x 10 -12 erg cm" 2 s" 1 , corrected for a dual absorber in the 
best fit. This differs only slightly from the unabsorbed flux of -F2-10 ~ 9x 10 _12 ergcm -2 s" 1 
we have obtained from our data in the present epoch. GOO use only one cold intrinsic 
absorber in their best fit, which yields an unabsorbed flux of F2-10 ~ 9 x 10~ 12 erg cm" 2 s" 1 , 
also consistent with our 2007 observation. 



The neutral hydrogen absorbing column in our best fit is less dense than that of W99 

0.04 
0.08 



or GOO. We find evidence for one intrinsic neutral absorber of ~ 1.08+qqo x 10 23 cm 



(with a partial-covering fraction of ~ 83+4%). In their best-fitting model, W99 employ a two- 
zone intrinsic absorption model with iV H (l) = 3.001^6 x 10 23 cm" 2 and iV H (2) = 4.9± 2 ;4 x 
10 22 cm -2 . The GOO best-fit model uses one intrinsic absorber of iV H = 21q; 5 x 10 23 cm -2 . 
Our model contains a component of partly- ionized absorption that W99 and GOO do not, 
however: the absori model representing this absorbing material has a column density of 
Nr = 1.37lo 78 x 10 22 cm -2 and an ionization parameter of £ = 68+34° erg cm -1 s" 1 . Even 
with this additional layer the W99 and GOO models contain higher absorbing columns than 
our 2007 model, however. 

The heavy absorption in W99, in particular, resulted in a power-law photon index split 
between intrinsic and scattered components, where a power- law photon index of V = 1.7 was 
assumed and 11% was scattered, while 70% of the intrinsic flux was absorbed by the larger 
column and 30% by the smaller column. We attempted such a neutral dual-absorber fit (see 
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Fig. 7. — The relative contributions of each model component in the three relevant fits to 
the NGC 1052 XIS and PIN data: laor (top, Model 1), pexrav+laor (middle, Model 2) and 
kdblur (ref lion) (bottom, Model 3). Each model presents its own strengths and weaknesses, as 
discussed in the text. Models are extended to 100 keV to illustrate their projected contributions at 
higher energies. Red represents the powerdaw (affected by Galactic and intrinsic absorption), green 
represents the thermal mekal component, Gaussian emission lines are in dark blue and reflection 
features such as the laor line (top), the pexrav component (middle) and the kdblur (ref lion) 
model (bottom) are in light blue. The laor line in Model 2 is in magenta. 
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Component 


Parameter Value 


Model 1 


Model 2 


Model 3 


mekal 


kT ( keV) 


64+°'° 4 




64 +0 - 04 

u - u ^-0.04 




64 + °- 04 

u - u ^-0.04 






K kT (ph cm _2 s _1 ) 


4.141};°° x 10- 


5 


4.14l°] 4 x 10- 


5 


4.14l°| 4 x 10" 


-5 




F kT ( erg cm' 2 s _1 ) 


9.0011;^ x 10- 


14 


9.4811;^ x 10" 


14 


8.781};^ x 10" 


-14 


absori 


r 

iV H (xl0 22 cm- 2 ) 
T(K) 

£ (erg cm" 1 s -1 ) 


1 O7+0.54 
J-- ' -0.78 

3 x 10 4 
68l^°° 




i.uu_ .05 
1 07+O.27 
l -° ' -0.23 

3 x 10 4 




-■-• du -0.05 
I 4^+0-63 
i -^ d -0.53 

3 X 10 4 
6711 




zpcf abs 


iV H (xl0 22 cm- 2 ) 
%cover 


10.82 + °- 7 4 
0.841°;°* 




10-821°;^ 
0.84l°;° 2 




10.83!°;!f 
0.85l°;° 2 




zpo 


r 


1-5018S 




1.50l°:°J 




i-50l°;°J 






Kr ( ph cm- 2 s _1 ) 


l-55±g;|2 x 10- 


3 


1.55l°; 2 5 x 10- 


3 


1.62l°;| x 10- 


•3 




F r (erg cm- 2 s -1 ) 


1.89l°; 24 x 10- 


11 


1.79l°; 24 x 10" 


11 


!-85la32 x 10" 


-11 


zgauss 


E(keV) 


fi 40+°- 01 

°-^ u -0.01 




6 40 +0 - 01 
°-^ u -o.oi 




6 40+ 001 
°-^ u -o.oi 






^narrow (Ph Cm" 2 S- 1 ) 


l-18±g;i| x 10" 


5 


i.i9±o:ig x 10- 


5 


l-26±g;i| x 10 ~ 


-5 




^narrow (erg Ctll" 2 ) 


l.OOlJjijj x 10 ~ 


13 


9.54l}| 4 x 10" 


14 


1.231°-^ x 10" 


-13 




EW narrow ( eV) 


1 1 1 + 18 

iii -18 




lH-il 




i2i±ig 




laor, kdblur 


C^emis 

r in (r g ) 
' (0) 


2-40l 2 ;^ 

19.931M 
72 +o 

' z -15 




2-40l 2 i 5 7 

19.76t£$g 

72l° 5 




o.gsio 1 ; ! 

on 1 £+25.36 
zu - id -20.15 
79+O 
' z -15 






Abroad (Ph Cm' 2 S- 1 ) 


l-69±5;g| x 10- 


5 


i-7i±8:£ x 10- 


5 








Abroad ( erg cm" 2 s -1 ) 


1.441°;™ x 10- 


13 


1.3818;™ x 10- 


13 








E^T XT ( ,.\ T\ 

E 14/broad (eV) 


1851^ 




187l 7 ° 6 








pexrav, reflion 


Fe / solar 

£ ( erg cm -1 s -1 ) 






1.0 
0.0 




1.0 

109l 2 09 






Rre& 






01 +0 
u - Ui -o.oi 




5.751° 75 x 10- 


-3 




K IcR (ph cm~ 2 s- 1 ) 






7.2ll° 21 x 10" 


8 


1.301};^ x 10" 


-7 




F rcfl (erg cm' 2 s- 1 ) 






4.041° 04 x 10" 


18 


4.24 x 10~ 13 






X 2 /" 


660/685 (0.96) 


660/683 (0.97) 


660/684 (0.96) 



Table 1: Comparison of our three best-fit models for the 0.5 — 35keV spectrum of NGC 1052. 
All components of Models 1-3 are modified by Galactic hydrogen absorption with Ah = 2.83 x 
10 20 cm . "K" denotes the normalization value of a given component. Flux values indicate 
absorbed flux from 0.5 — 35 keV. Redshifts for the model components are held constant at the 
cosmological value for NGC 1052: z = 0.0049. Abundances not listed are held constant at their 
solar values. The inclination angle of the accretion disk used in the laor and kdblur models was 
constrained to fall within the radio observation uncertainties of i = 57 — 72°. All errors listed are 
at 90% confidence for one interesting parameter. 



Waveband (keV) 



Model 1 Flux ( erg cm 2 s 1 ) 



3.42 x lO" 13 
3.42 x 1(T 12 



Model 2 Flux ( erg cm" 2 s" 1 ) 



3.42 x 10" 13 
3.42 x 1(T 12 



Model 3 Flux ( erg cm" 2 s" 1 ) 



0.5-2.0 



3.68 x 10~ 13 
5.31 x 10^ 12 



5.53 x 10" 12 
8.96 x 10~ 12 



2.0 - 10.0 



5.25 x 10~ 12 
8.96 x 10~ 12 



5.33 x 10~ 12 
9.20 x HT 12 



10.0-60.0 



3.64 x 10~ n 
3.68 x 10- 11 



3.64 x 10~ n 
3.68 x 10- 11 



3.79 x 10~ n 
3.82 x 10- 11 



to 

03 



Table 2: Extrapolated model fluxes for the 2007 Suzaku observation of NGC 1052 in three wavebands. Absorbed fluxes are on 
the top line, unabsorbed fluxes are on the line below. For the softest energies the flux from XIS 1 was used, while the XIS 0+3 
data was used for the 2 — lOkeV band. PIN data were used for the flux above lOkeV. 
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§4.1114.21) and found that the parameters of the second absorber could not be constrained, 
and that a single intrinsic neutral absorber and an ionized absorber along with Galactic 
photoabsorption modifying a simple, albeit fairly hard power-law (T ~ 1.5) provided an 
excellent fit to the continuum. 

W99 and GOO also found evidence for a thermal mekal component to explain the soft 
excess noted below 2 keV, which was best modeled with kT = 0.53^ keV and a 0. 1-2.0 keV 
flux of .F m ekai = 5.8+Jjj x 10 -14 erg cm -2 s _1 in W99. GOO modeled this component with a 
similar model of much higher temperature of kT > 5keV, th ough these authors d o not 



report its flux. This thermal component was also detected by iKadler et al.l (j2004al ) with 



Chandra at kT = 0.41+°;°?, though over a flux range of 0.2 — 8.0 keV and with limited 
signal-to-noise. We found a similar thermal component in the Suzaku data. Though the XIS 
detectors lose effective area rapidly below 0.3 keV, extrapolating Model 1 down to 0.1 keV 
yielded kT = 0.64 ± 0.04 keV and F mekal = 9.38+Jjg x 10- 14 ergcm- 2 s" 1 over 0.1 - 2.0 keV. 
Our results showed a slight increase in temperature and flux for the mekal component but 
were consistent with the 1996 results within error bars. 

In the favored dual-absorber ASCA model of W99, a broad Fe Ka line component was 
not required in the fit and hence these authors found no definitive evidence for reflection 
from an optically-thick accretion disk in their data. They did require a narrow Fe Ka line, 
however: the centroid energy of this component was held constant at E — 6.37 keV and 
had an equivalent width of EW = 270 ± 120 eV. GOO did detect an Fe K line with some 
apparent broadening, though this could be accounted for in the BeppoSAX data with a 
mildly ionized Fe Ka line with centroid energy E = 6.48+°'2okeV and an equivalent width 
of EW = 230 ± 170 eV. In our Model 1 for the Suzaku data, we found that both broad and 
narrow components were necessary in order to eliminate the residuals between 4 — 7keV. Our 
narrow line, with a centroid energy of 6.4 ± 0.01 keV had an equivalent width of 111 ± 18 eV. 
The presence of a broad component to the Fe Ka line likely rendered this narrow feature 
significantly weaker than its 1996 counterpart. 

Our broad Fe Ka line, held constant at a centroid energy of E = 6.4 keV, was best fit 
with a laor emission feature originating from the accretion disk around a rapidly rotating 
black hole, though the inner radius of emission within the disk was only mildly constrained 
to ri n < 45 r g , as discussed in §4.21 This feature was stronger than the narrow Fe Ka line: 
its equivalent width was 185±75 2 eV. Replacing the laor line with a broad Gaussian, we 
found that a = 1.04 ± 0.26 keV, which translates to a lower-limit FWHM of v ~ 0.37c and 
implies an origin for this feature very close to the black hole. Though the fit of the laor line 
itself is suggestive of a broad, fluorescent Fe Ka origin in the inner accretion disk, the lack 
of reflected continuum flux above lOkeV argues against this interpretation, as we discuss 
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below and in the following Section. 

W99 considered a reflection-based href 1 model to fit the ASCA spectrum, but were 
unable to obtain meaningful constraints on the reflection parameters without making several 
assumptions. These included fixing the value of the power-law photon index of the central X- 
ray source to be V = 1.7 and also fixing the inclination angle of the disk to i = 60°. The latter 
was consistent with our own angle constraints based on the VLBI observations of NGC 1052, 
but our Suzaku observation indicated a somewhat harder photon index (see Tabled]). For an 
assumed reflection fraction of -R re fl = 1 and the observed ASCA spectrum, W99 found that 
only a small fraction of the direct emission from the X-ray source was visible, leading these 
authors to infer that a thick absorbing column in the nucleus (Nn > 3 x 10 23 cm -2 ) blocked 
the X-ray source from view, such that reflected emission dominated the spectrum. However, 
the equivalent width of the Fe Ka line, EW ~ 40 eV, was over an order of magnitude 
smaller than the same line seen in other heavily obscured AGN. To account for this line with 
a reflection model would require an implausibly small iron abundance of Zp e /solar < 0.05. 
This reflection model was therefore argued against by W99. Similarly, GOO also ruled out 
evidence for Compton reflection in the NGC 1052 BeppoSAX spectrum due to a statistically 
unacceptable fit for the pexriv model. These authors were only able to mildly constrain the 
reflection fraction to a hard upper limit of _R re fi < 0.6. 

The 2007 Suzaku data concur with this assessment of a lack of reflection: though Model 2 
produced a comparable global goodness-of-fit to Model 1, the reflection fraction was found 
to be quite small: _R rc fi < 0.01. Attempting to model any reflection present with Model 3, we 
were unable to constrain the emissivity index of the accretion disk, though we did achieve 
similar constraints for the equivalent width of the narrow iron line and the photon index 
of the power-law component. Based on this spectral fitting, we conclude that in 2007, 
NGC 1052 does show evidence for a broadened iron line; this line is not particularly broad 
or strong, however, and is not accompanied by the expected Compton hump above 10 keV. 
This may indicate that the line arises outside of the innermost portions of the disk and/or 
that reflection simply does not play a large role in creating the overall X-ray spectrum. 

In order to make a more detailed comparison between the 1996 ASCA data and the 2007 
Suzaku data for NGC 1052, we have examined the 0.5 — 10 keV W99 spectrum in combination 
with our Model 1. Data from the SIS and GIS 2 were used, as these were the instruments 
with the highest count rates and the closest on-axis telescope pointing during the ASCA 
observation. We show a plot of the ratio of the ASCA data to our Suzaku best-fit Model 1 
in Fig. E 

Changes between the two epochs are evident at energies below ~ 3keV and above 
~ 7keV, likely indicating that the continuum and absorbing structures have undergone 
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Fig. 8.— Ratio plot of the 1996 ASCA spectrum of NGC 1052 to the 2007 Suzaku best-fit Model 1. 
SIS and GIS 2 data are shown with black and red crosses, respectively. The green line denotes a 
theoretical perfect fit. Note the curvature of the residuals below ~ 3keV. These residuals indicate 
that changes have occurred in the absorption structure surrounding the AGN between epochs. 
x 2 / v = 255/141 (1.81). No refit has taken place. 
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physical changes during the intervening years. Refitting the model, we found that x 2 l v — 
120/141 (0.85) as opposed to x 2 / u — 255/141 (1.81) before refitting. Reduced chi-squared 
values less than unity often indicate that the model has more parameters than are necessary 
to adequately fit the data, and given the lower spectral resolution and throughput of the 
ASCA data, this is likely the case here. Nonetheless, we can use this refit to verify which 
model components have undergone significant changes between 1996 and 2007. 

The thermal component had a temperature of kT = 0.50 ± 0.14keV, consistent with 
the results of W99 and consistent in its upper limit with Model 1. The normalizations of 
the thermal component were also consistent within error bars. The power-law photon index 
from Model 1 was T = 1.50t^, vs. T = 1.28 ± 0.91 from our refit of the ASCA data, also 
consistent within errors. W99 had maintained this component at a fixed value of T = 1.7 
for their best-fitting models, but our refitted value was consistent with the W99 fits for 
their models 2-4. Though W99 employed a power-law component split between intrinsic and 
scattered emission while we determined that only intrinsic emission was necessary for the 
Suzaku data, the total normalization of all the W99 power-law components was consistent 
with our intrinsic power-law normalization within errors. Our refitted ASCA value (assuming 
intrinsic emission only) was also consistent with the W99 intrinsic normalization within 
errors. Refitting the absorption parameters yielded an ionized absorbing column of Nn < 
5.64 x 10 21 cm~ 2 and an ionization parameter of £ < 97 erg cms" 1 . While these values 
were not well-constrained, they were at least marginally consistent with those of Model 1. 
However, while the partial-covering absorber column was constrained to Nn < 15.19 x 
10 22 cm -2 , consistent with Model 1, the covering fraction could not be constrained. These 
changes suggest that the absorbing structure differs physically in 2007 as compared with 
1996. 

The iron line complex, by contrast, showed little variation between the two epochs. 
Our refit to the ASCA data yielded a narrow Fe Ka component of normalization K = 
1.19io/78 x 10~ 5 phcm" 2 s -1 , consistent with our Suzaku Model 1 value and the W99 fit, and 
strongly indicating that this component, if variable, is only variable on long timescales and 
must therefore be emitted from material far from the black hole. W99 reported that no broad 
Fe Ka line component was robustly present in 1996, contrary to our findings in 2007. Our 
refit, including a laor component, does not adequately constrain the parameters of the broad 
line in 1996, confirming the null result of W99. This could mean either that the broad line 
component was truly absent in 1996, or perhaps that it was present at a lower normalization 
and effectively drowned out by the continuum emission. Given the rapid decrease in effective 
area above ~ 5keV and the spectral resolution of ASCA, however, it is very likely that the 
broad iron line component simply could not have been detected by this observatory, thus 
rendering moot any comparison between a broad Fe Ka line in 1996 and 2007. 
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6. Discussion 



We have performed detailed spectral fits to the XIS and HXD/PIN data from the 101 ks, 
2007 Suzaku observation of the LINER galaxy NGC 1052. Three different spectral models 
were employed: (1) an absorbed power-law continuum with contributions from a soft thermal 
component along with narrow and broad Fe Ka lines; (2) the same model, but with an added 
component of reflected continuum emission from a neutral accretion disk (pexrav); (3) the 
same base model, again, but with a more advanced, self-consistent reflection model which 
includes fluorescent emission lines from the disk and takes ionization into account fref lion; 



sec 



Ross fe Fabianl (120051 ) for details). In Model 3 we applied relativistic smearing using the 



kdblur model. Because fluorescent lines are included in ref lion, there is no need for an 
additional broad iron line component as in Model 1. We summarize our findings below: 



• The three models are quite similar in their statistical goodness-of-fit, though some of 
the disk reflection parameters in Models 2-3 prove especially difficult to constrain. 

• The continuum is well described by a component with a power-law spectrum of photon 
index T ~ 1.5, thermal emission with kT ~ 0.64 keV, a neutral intrinsic absorbing 
column of Ah ~ 10 23 cm -2 and a covering fraction of ~ 84%, and an ionized intrinsic 
absorber with A^h ~ 1-37 x 10 22 cm -2 and £ ~ 68ergcm~ 1 s -1 . These findings are 
relatively consistent with those of the 1996 ASCA observation of W99, though our 
photon index is harder and our absorbing column lower in density in 2007. 

• As in the ASCA data, a narrow Fe Ka line is present: EW ~ llleV. 

• Strong evidence exists for a broadened Fe Ka emission component (EW ~ 185 eV). 
Though the best statistical fit is achieved with a laor line profile, however, the rather 
broad constraint on the inner radius of r in < 45 r g renders it impossible to distinguish 
between a spinning and non-spinning black hole. 

• The PIN instrument did detect hard X-ray emission above lOkeV from NGC 1052. 
This emission is well fit by the power-law component of the 2 — lOkeV continuum 
and does not show any evidence for reflected emission via the Compton hump that is 
commonly associated with other AGN harboring broad iron lines, and expected based 
on the presence of the broad iron line if the two features arise from reflection onto a 
standard accretion disk. Neither of our reflection models found any statistical evidence 
for a strong reflection component. 

• Though VLBI observations of the jet in NGC 1052 have constrained the orientation 
angle of the disk to be between i = 57 — 72°, we find that removing this constraint 
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in our fitting yields a best-fit inclination angle of i = 45 ± 5°, though this does not 
significantly improve the global goodness-of-fit. Nonetheless, this result does mark a 
~ 2o" deviation from the radio results. 

These findings motivate three main questions: (1) what is the physical origin of the 
broad iron line in our 2007 observation; (2) how is it that we detect a broad iron line yet 
no evidence for disk reflection above lOkeV, and (3) why is the laor line inclination, when 
kept as a free parameter in the fit, not more consistent with the disk inclination constraints 
established in radio observations of the inner jet in NGC 1052? 

As it turns out, the answers to these questions may be closely related. Relativistically 
broad iron lines are expected to be associated with pronounced Compton humps. Comp- 
tonized reflection originates via irradiation of the optically thick inner accretion disk by 
the hard X-ray source, perhaps in some sort of lamp-post geometry, e.g., the base of a 



jet ( jMiniutti et al.l 120071 ). This same process will also produce fluorescent line emission 
from the disk, most notably the Fe Ka line, which will experience significant broaden- 
i ng and skewness by virt ue of its origin from the spacetime proximal to the black hole 



([Reynolds fe Nowata 120031 ) . The absence of a Compton hump, however, casts doubt on the 
presence/contribution of reflection from the disk to the overall X-ray spectrum. From our 
pexrav model fit (Model 3), we can estimate the equivalent width of the iron line we should 
expect if both the Compton hump and the broad line originate from disk reflection. Using 
the ratio of the normalization of the power-law to the upper limit of the reflected compo- 
nent in our fit, we infer that the iron line equivalent width should be EW ~ 80 eV if the 
two features are indeed produced by the same physical process. This low value is highly 
inconsistent with the observed strength of the broad iron feature, however. 

Relativistic disk smearing is not the only process by which this spectral line can be 
broadened. In the following paragraphs we consider a number of other physical scenarios, 
including: transition to an advection-dominated accretion flow (ADAF) in the inner disk, 
resonant scattering from the putative torus region, origin of the line in an outflow located 
in the broad line region (BLR), and Fe Ka produced by Cerenkov line-like radiation rather 
than fluorescence. 

One explanation for the lack of an observed reflected continuum could be that the accre 



tion d isk transitions to a radiatively-inefficient state within some radius, e.g.. iNarayan fe Yi 



(If994j ). In this ADAF picture, which has often been invoked as a potential explanation for 
low-luminosity sources (especially those with observed jet activity), the disk traps radiation 
rather than releasing it, puffing up to become optically thin and geometrically thick as its 
temperature and ionization state rises. Because of these properties, the gas in this region 
would not be expected to produce much, if any, reflection, either in the form of a Compton 
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hump or discrete emission lines. If a contribution to the broad iron line did somehow arise 
from this region, we would expect it to be highly ionized, yet this is not seen in the data: 
our broad and narrow iron line components both have energies robustly less than 6.41 keV, 
suggesting that they arise from species at or below Fe xvill. We therefore conclude that an 
ADAF, if present, is unlikely to be a significant contributor to the broadening of the Fe Ka 
feature. 

The lack of a Compton hump at the expected strength and the mismatch between the 
unconstrained laor disk inclination and the VLBI value causes us to rule out a relativistic 
inner disk origin for the broad iron line. We also eliminate a possible origin in an ADAF 
due to ionization constraints and optical depth issues. Investigating other possible sources 
of iron emission in the central AGN leads us to consider the role of resonant scattering 
in the torus region, which is thought to contribute heavily to the flux of the narrow Fe 
Ka line. However, both the large equivalent width and the large FWHM velocity of the 
broad line component we observe (both substantially greater t han the v ~ 2500 km s - 1 



and EW < 100 eV typically associated with the narrow core, e.g., lYaqoob fe Padmanabhan 
(120041 )) render this explanation implausible. 



An AGN outflow originatin g in the BLR is exp lored as a possible explanation for the 
similar spectrum of NGC 7213 (IBianchi et al.ll2008l ). Here a broad iron line with no corre- 
sponding Compton hump is also observed, and these authors fi nd that the eq u ivalen t width 



of the Fe Ka feature is consistent with predictions made by lYaqoob et al 



(|200lh for an 



origin in the BLR or perhaps a Compton-thin torus. Further, measurements of the velocity 
broadening of the optical Ha line concur with those measured from the broad Fe Ka line 
(v = 2500 km s -1 ), suggesting a similar origin for the emission. Using the same diagnos- 
tic, we estimate that we should expect a broad line EW ~ 75 — 100 eV if it originates in 
an outflow from this region; this underestimates our observed equivalent width by nearly a 
factor of ten using a simple Gaussian line to represent the broad feature, and by roughly a 
factor of two if the better-fitting (though, in this case, unphysical) laor line is used. Also, 
our estimated lower-limit FWHM velocity of the line (v ~ 0.37c) is more th an an order 



of ma gnitude greater than the polarized optical Ha FWHM measurement of iBarth et al. 
( 119991 ) (v ~ 5000 km s -1 ). This is a strong indication that the broad Fe Ka line we observe 



originates closer to the black hole than th e BLR. Moreover, variability of the b road iron line 
on years- long timescales has been noted (IKadler et al.l 12004a : iRos et al.l 120071 ); this type of 
variation would not be expected from the typical symmetric BLR, as changes in flux due to 
passing clouds in this region should theoretically average out over time. 

Finally, we consider a more exotic possibility: perhaps the broad iron line i s prod uced 
by Cerenkov line-like radiation in the centralmost regions of the AGN. I You et al.l (120031 ) put 
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forward the intriguing notion that in dense regions where the refractive index of the material 
is large, relativistic electrons impinging upon this material can produce Cerenkov radiation 
in a narrow wavelength range very close to the intrinsic atomic wavelength of the material. 
The combination of absorption and emission causes the final emission feature to be redshifted 
and skewed, making it appear very like a relativistic line produced by fluorescence in the 
inner accretion disk. Because electrons are producing the emission rather than photons, no 
reflected continuum is expected. The cool, dense, iron-rich gas in the disk of NGC 1052 
could provide the medium for this process, while relativistic electrons in the corona or near 
the base of the jet would make for an ideal bombardment population. Furthermore, this 
physical picture makes for a much more consistent match with our estimate of the distance 
of the broad line region of origin from the hard X-ray source: our lower-limit FWHM velocity 
for the broad Gaussian fit to the line yields v ~ 0.37c, which corresponds to an upper-limit 
distance of d ~ 8 r g . This is well inside the BLR. 

Given our constraint on the distance of the broad iron line emitter from the hard X-ray 
source, iron fluorescence from the base of the jet itself is also a potential candidate for the 
emission mechanism. However, considering the orientation of the jets observed in NGC 1052 
to our line of sight, it is unclear whether any fluorescent emission from the jet base could 
escape the system without being absorbed. Correlated radio flux measurements from the 
base of the jet and X-ray spectra of this source must be examined over several epochs, 
simultaneously whenever possible, in order to assess th e likelihood of th is scenario as well 
as that of the Cerenkov line-like radiation proposed by I You et al.l (120031 ) . Only through a 
coordinated study of the jet and the accretion flow in NGC 1052 can we hope to understand 
the connection between these two vital physical processes in AGN. 



7. Conclusions 

Our 2007 Suzaku spectrum of the LINER galaxy NGC 1052 is consistent with X-ray 
spectra of this source from previous epochs, with a fairly flat power-law continuum that is 
heavily absorbed. Intrinsic neutral and ionized columns are detected, along with evidence of 
Galactic photoabsorption. A thermal component is also detected, likely due to the interaction 
of the jets with the surrounding ISM. Both narrow and broad iron lines are observed, though 
interestingly there is little to no reflection seen in the spectrum above lOkeV: two different 
models used to characterize this feature both require i? re fl < 0.01. 

We are thus faced with a complex scenario for NGC 1052 in which the broad iron line 
is not consistent with the lack of observed continuum reflection from the disk. While it is 
possible that the reflection has been drowned out in some way by the power-law continuum, 
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we must also consider the physical implications if reflection is indeed absent above 10 keV. It 
could be that we are witnessing the transition of the accretion flow to a radiatively-ineflicient 
state at some critical radius in the disk, which could largely eliminate reflection features from 
the spectrum. A broad iron line is still produced, however, and the velocity width of the 
line is consistent with an origin close to the black hole. If the inner disk is an ADAF, it is 
possible that the broad Fe Ka line is emitted from optically-thick material near the base of 
the jet, or per haps from a yet m ore exotic mechanism such as the Cerenkov-like line radiation 
postulated by You et al. ( 20031 ). Coordinated observations of both the inner jet(s) (radio) 
and the inner accretion flow (X-ray) are needed in order to solve this puzzle. 



Acknowledgments 

LB and MK thank the NASA Postdoctoral Program, administered by ORAU, for 
their support at NASA's GSFC. KW gratefully acknowledges support from NASA grant 
NNX08AC226, MA from NASA grant NNX08AC22G. YYK is a Research Fellow of the 
Alexander von Humboldt Foundation. LB and KW also appreciate excellent advice from 
Tahir Yaqoob, Chris Reynolds, Julian Krolik, Demos Kazanas and Tim Kallman. KW would 
like to dedicate this work to Professor Andrew Wilson, who prompted her interest in LINERs 
and NGC 1052 in the 1990s. We gratefully acknowledge the helpful comments received from 
our anonymous referee, which have improved this work. 



REFERENCES 

Arnaud, K. A. 1996, in Astronomical Society of the Pacific Conference Series, Vol. 101, 
Astronomical Data Analysis Software and Systems V, ed. G. H. Jacoby & J. Barnes, 
17-+ 

Barth, A. J., Filippenko, A. V., & Moran, E. C. 1999, Astrophys. J. Lett, 515, L61 

Bianchi, S., La Franca, F., Matt, G., Guainazzi, M., Jimenez Bailon, E., Longinotti, A. L., 
Nicastro, F., & Pentericci, L. 2008, Mon. Not. R. astr. Soc, L74+ 

Blandford, R. D. & Znajek, R. L. 1977, Mon. Not. R. astr. Soc, 179, 433 

Braatz, J. A., Wilson, A. S., & Henkel, C. 1994, Astrophys. J. Lett, 437, L99 

Claussen, M. J., Diamond, P. J., Braatz, J. A., Wilson, A. S., & Henkel, C. 1998, Astro- 
phys. J. Lett, 500, L129+ 



-33- 



Fabian, A. C, Rees, M. J., Stella, L., & White, N. E. 1989, Mon. Not. R. astr. Soc, 238, 
729 

Gabel, J. R., Bruhweiler, F. C, Crenshaw, D. M., Kraemer, S. B., & Miskey, C. L. 2000, 
Astrophys. J., 532, 883 

Gierlihski, M., Done, C, & Page, K. 2008, MNRAS, 388, 753 

Guainazzi, M. & Antonelli, L. A. 1999, MNRAS, 304, L15 

Guainazzi, M., Oosterbroek, T., Antonelli, L. A., & Matt, G. 2000, Astr. Astrophys., 364, 
L80 

Hawley, J. F. & Krolik, J. H. 2006, ApJ, 641, 103 
Heckman, T. M. 1980, Astr. Astrophys., 87, 152 

Ho, L. G, Filippenko, A. V., & Sargent, W. L. W. 1997, Astrophys. J., 487, 568 

Kadler, M., Kerp, J., Ros, E., Falcke, H., Pogge, R. W., & Zensus, J. A. 2004a, Astr. Astro- 
phys., 420, 467 

Kadler, M., Ros, E., Lobanov, A. P., Falcke, H., & Zensus, J. A. 2004b, Astr. Astrophys., 
426, 481 

Kadler, M., Ros, E., Weaver, K., Kerp, J., & Zensus, J. A. 2004c, in Bulletin of the American 
Astronomical Society, Vol. 36, Bulletin of the American Astronomical Society, 823 — h 

Kalberla, P. M. W., Burton, W. B., Hartmann, D., Arnal, E. M., Bajaja, E., Morras, R., & 
P6ppel, W. G. L. 2005, Astr. Astrophys., 440, 775 

Kataoka, J., Reeves, J. N., Iwasawa, K., Markowitz, A. G., Mushotzky, R. F., Arimoto, M., 
Takahashi, T., Tsubuku, Y., Ushio, M., Watanabe, S., Gallo, L. G, Madejski, G. M., 
Terashima, Y., Isobe, N., Tashiro, M. S., & Kohmura, T. 2007, Publ. astr. Soc. Japan, 
59, 279 

Kinkhabwala, A. A. 2003, PhD thesis, AA(COLUMBIA UNIVERSITY) 
Knapp, G. R., Faber, S. M., & Gallagher, J. S. 1978, Astr. J., 83, 139 
Laor, A. 1991, Astrophys. J., 376, 90 

Magdziarz, P. & Zdziarski, A. A. 1995, Mon. Not. R. astr. Soc, 273, 837 
Markowitz, A. & Edelson, R. 2004, Astrophys. J., 617, 939 



-34- 



Marscher, A. P. 2006, Astron. Nachr., 327, 217 

Marscher, A. P., Jorstad, S. G., Gomez, J.-L., Aller, M. F., Terasranta, H., Lister, M. L., & 
Stirling, A. M. 2002, Nature, 417, 625 

Meier, D. L., Koide, S., & Uchida, Y. 2001, Science, 291, 84 

Mewe, R., Gronenschild, E. H. B. M., & van den Oord, G. H. J. 1985, Astr. Astrophys. Suppi, 
62, 197 

Miller, J. M., Homan, J., Steeghs, D., Rupen, M., Hunstead, R. W., Wijnands, R., Charles, 
P. A., & Fabian, A. C. 2006, ApJ, 653, 525 

Miniutti, G., Fabian, A. C, Anabuki, N., Crummy, J., Fukazawa, Y., Gallo, L., Haba, 
Y., Hayashida, K., Holt, S., Kunieda, H., Larsson, J., Markowitz, A., Matsumoto, 
C, Ohno, M., Reeves, J. N., Takahashi, T., Tanaka, Y., Terashima, Y., Torii, 
K., Ueda, Y., Ushio, M., Watanabe, S., Yamauchi, M., & Yaqoob, T. 2007, 
Publ. astr. Soc. Japan, 59, 315 

Narayan, R. & Yi, I. 1994, Astrophys. J. Lett, 428, L13 

Nelson, C. H. & Whittle, M. 1995, ApJS, 99, 67 

Papadakis, I. E., Petrucci, P. O., Maraschi, L., McHardy, I. M., Uttley, P., & Haardt, F. 
2002, Astrophys. J., 573, 92 

Rees, M. J., Begelman, M. C, Blandford, R. D., & Phinney, E. S. 1982, Nature, 295, 17 

Reynolds, C. S. & Nowak, M. A. 2003, Phys. Reports, 377, 389 

Ros, E., Kadler, M., Kaufmann, S., Kovalev, Y. Y., Tueller, J., & Weaver, K. A. 2007, in 
Highlights of Spanish Astrophysics IV, ed. F. Figueras, J. M. Girart, M. Hernanz, & 
C. Jordi, 165-+ 

Ross, R. R. & Fabian, A. C. 2005, Mori. Not. R. astr. Soc., 358, 211 
Shakura, N. I. & Syunyaev, R. A. 1973, Astr. Astrophys., 24, 337 

Tremaine, S., Gebhardt, K., Bender, R., Bower, G., Dressier, A., Faber, S. M., Filippenko, 
A. V., Green, R., Grillmair, C, Ho, L. C, Kormendy, J., Lauer, T. R., Magorrian, 
J., Pinkney, J., & Richstone, D. 2002, ApJ, 574, 740 

Weaver, K. A., Wilson, A. S., Henkel, C, & Braatz, J. A. 1999, Astrophys. J., 520, 130 



-35 - 

Woo, J.-H. & Urry, C. M. 2002, ApJ, 579, 530 

Yaqoob, T., George, I. M., Nandra, K., Turner, T. J., Serlemitsos, P. J., & Mushotzky, R. F. 
2001, Astrophys. J., 546, 759 

Yaqoob, T. & Padmanabhan, U. 2004, Astrophys. J., 604, 63 

You, J. H., Liu, D. B., Chen, W. P., Chen, L., & Zhang, S. N. 2003, Astrophys. J., 599, 164 

Zdziarski, A. A., Johnson, W. N., Done, C, Smith, D., & McNaron-Brown, K. 1995, Astro- 
phys. J. Lett, 438, L63 

Zhang, Y. H., Treves, A., Maraschi, L., Bai, J. M., & Liu, F. K. 2006, Astrophys. J., 637, 
699 



This preprint was prepared with the AAS IATgX macros v5.2. 



